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Abstract. The purpose of this note is to establish a Central Limit Theorem 
for the number of eigenvalues of a Wigner matrix in an interval. The proof 
relies on the correct aymptotics of the variance of the eigenvalue counting 
function of GUE matrices due to Gustavsson, and its extension to large families 
of Wigner matrices by means of the Tao and Vu Four Moment Theorem and 
recent localization results by Erdos, Yau and Yin. 



1. Introduction 

This note is concerned with the asymptotic behavior of the eigenvalue counting 
function, that is the number Nj of eigenvalues falling in an interval /, of families 
of Wigner matrices, when the size of the matrix goes to infinity. Wigner matrices 
are random Hermitian matrices M„ of size n such that, for i < j, the real and 
imaginary parts of {Mn)ij are iid, with mean and variance ^, {Mn)ii are iid with 
mean and variance 1. An important example of Wigner matrices is the case where 
the entries are Gaussian, giving rise to the so-called Gaussian Unitary Ensemble 
(GUE). GUE matrices will be denoted by M'^. In this case, the joint law of the 
eigenvalues is known, allowing for complete descriptions of their limiting behavior 
both in the global and local regimes (cf. for example [1]). 

Denote by Ai,...,A„ the real eigenvalues of the normalized Wigner matrix 
Wn = ^M„. The classical Wi gner theorem states that the empirical distribution 

i X]J=i '^Aj on the eigenvalues of Wn converges weakly almost surely as n — )• oo to 
the semicircle law dpsdx) = — x'^l[-2,2]{x)dx. Consequently, for any interval 

/CM, 



1 1 
^/(W„) = - 1{A,.G/} Psc{I) almost surely. 



At the fluctuation level, it is known, due to the particular determinantal structure 
of the GUE, that 

Theorem 1 (Costin-Lebowitz g], Soshnikov [7] (see [J)). Ltt be a GUE ma- 
trix. Set = -i=M'„. Let In he an interval in M. // Var(iV/„ (IV,' )) oo, 

V ^ n— >oo 
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then 

(1) N,AW:^)-nNrAW^)] ^ 

^Var(Ar/„(H/;j) 

in distribution. 



In 2005, Gustavsson [5] was able to fully describe for GUE matrices the as- 
ymptotic behavior of the variance of the counting function Ni{W^) for intervals 
/ = [y, +oo) with y E (—2, 2) strictly in the bulk of the semicircle law. He estab- 
lished that 

(2) E[iV,«)] = npsciI) + 0[^) and Var(Ar,(t^;;)) = (^+o(l)) logn. 
In particular therefore, if / = [y, +oo) with y e (—2, 2), 

(3) N,iw:^)-miiwD] ^ 

as well as 

(4) mWl^)^npUl) ^ ^(0,1) 

^2¥^logn 

(which we call below the CLT with numerics). 

The purpose of this note is to extend these conclusions to non-Gaussian Wigner 
matrices. The class of Wigner matrices covered by our results is described by the 
following condition. Say that Af„ satisfies condition (CO) if the real part ^ and the 
imaginary part ^ of (M„)ij are independent and have an exponential decay: there 
are two constants C and C such that 

P(|ei > t^) < e-' and > t^) < e~\ 

for aU t>C'. 

Say that two complex random variables ^ and ^' match to order k if 

E [Re(e)'"Im(C)'] = E [Re(C')™Ini(e')'] 
for all TO, Z > such that ni + I < k. 

The following theorem is the main result of this note. 

Theorem 2. Let A/„ be a random Hermitian matrix whose entries satisfy condition 
(CO) and match the corresponding entries of GUE up to order 4. Set Wn — "^-^'-^n- 
Then, for any y G (—2,2) and I{y) = [y, +oo), setting y„ iV/(j,)(M^„), we have 

E[r„] = npsc{I{y)) + o(l) and Var(r„) = + o(l)) logn, 

and the sequence (Yn) satisfies the CLT 

~ ^^^"^ AA(0,1). 
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The theorem is estabhshed in the next two sections. In a first step, relying on 
Gustavsson's resuhs and its extension to Wigner matrices by Tao and Vu [9], we 
estabhsh that satisfies the CLT with numerics (jH). In a second step, we use 
recent results of Erdos, Yau and Yin 3 on the localization of eigenvalues in order 
to prove that E[y„] and Var(l^) are close to those of M^^ (GUE) and therefore 
satisfy 



On the basis of the CLT with numerics, Gustavsson [5] described the Gaussian 
behavior of eigenvalues in the bulk of the semicircle law in the form of 

(5) pWA.(^-)_t(Vn) ^ 

n 

in distribution, where t{x) e [—2, 2] is defined for x e [0, 1] by 



More informally, \r{Wl^) « t(i/n)+7V(0, (4_t^( ■/^"2)„2 )■ This is achieved by the tight 
relation between eigenvalues and the counting function expressed by the elementary 
equivalence, for I{y) = [y, +oo), y G M, 



The result ([5]) was extended in [9^ to large families of Wigner matrices satisfying 
condition (CO) by means of the Four Moment Theorem (see [9] and [10]). Now 
using the reverse strategy based on ([S]), ([S]) may be shown to imply back the CLT 
with numerics (|3]) for Wigner matrices whose entries match those of the GUE up 
to order 4. We provide some details in this regard relying on the following simple 
consequence of the Four Moment Theorem. 

Proposition 3. Let M„ and Af " be two random matrices satisfying condition (CO) 
such that their entries match up to order 4. There exists c > such that, if n is 
large enough, for any y £ (—2,2) and any (i,j) G {1, . . . if I{y) — [y, +cxd), 



2. CLT WITH NUMERICS AND EIGENVALUES IN THE BULK 




(6) 



Nji^y-fCWn) < n — i if and only if A.^ < y. 



P(A.e/(y))-P(A:'e/(2/))|<n 



and 



P(A, G I{y) A A, e I{y)) - P(A'/ e I{y) A A^' e /(y))| < n'^ 



As announced, we would like to show that the behavior of eigenvalues in the bulk 
(O extended to Wigner matrices leads to the CLT with numerics for such matrices, 
namely. 
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in distribution for Wigner matrices Wn satisfying (CO). To prove this, observe that 
for every x S M. 

' Nny-,{Wn)-npsc{I{y)) ^ \ 

— — , < a:^ I = P[Nj(^y)[W„) <n-i„) 



where i„ = npsc{{-oo, y]) - ^ logn. Then, by 

'Ni(y){W r,)-npsc{ I{y)) 

2^ logn 



V / 1 / 



V 2 S -^"J' 



where x„ = ^/ ^"^t'W") ^'^"^ Now ^ ^ p,c((-oo,y]) G (0,1). Furthermore, 
a;„ — ^ X since 



i(in/n) = t( p5c((-oo,j/]) - -W-i^logn 

n V Stt^ 



t (p,,((-(x), 2/])) - t' (p.e((-oo, y]))^y^logn + 



2 Vlogn , ,^Vlogn^ 
Hence = x^/ ^j^^^ + o(l), from which Xn — x. 



4 — j/^ n V n 



Applying ([5]) (extended to Wigner matrices), we obtain that 



V V 2 \/iQgyt / , - - 

n 

where X J\f{0, 1), implying d?]). 



3. Estimating the mean and the variance of 



To reach the CLT of Theorem[2]from the CLT with numerics (l7|) , it is necessary to 
suitably control the expectation and variance E[Yn] and Var(yri) of the eigenvalue 
counting function, and to show that their behaviors are identical to the ones for 
GUE matrices. The direct use of the Four Moment Theorem is unfortunately 
not enough to this purpose since it only provides proximity on a small number of 
eigenvalues. But recent results of Erdos, Yau and Yin [3], presented in the following 
statement, describe strong localization of the eigenvalues of Wigner matrices which 
provides the additional step necessary to complete the argument. 

Theorem 4. Let A/„ be a random Hermitian matrix whose entries satisfy condition 
(CO). There is a constant C > such that, for any i £ {1, . . . , n}, 

P{\X,-t{i/n)\ > (logn)^'°si°s"min(i,n-i + l)-i/3^-2/3) < ^-3^ 
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Note that if ne < i < (1 — e)n for some small e > 0, then min(i, n — i + 1) > ne 
so that 

(8) P(|A, - t{i/n)\ > 7i-ie-i/3(iog„)Ciogiogn^ < ^-3^ 

The next lemma presents the main conclusion on the expectation and variance 
of the eigenvalue counting function, extending Gustavsson's conclusion ([2]) for the 
GUE to Wigner matrices of the class (CO). Once this lemma is established, Theorem 
El will follow. 

Lemma 5. Set Wn = -^Mn, I{y) = [y,+oo) where y e (—2,2), and Yn = 
N^y^iWn). Then 

E[y„] = npsc{I{y)) + o(l) and Var(r„) = + o(l)) logn. 

Proof. By Gustavsson's results ^ therefore, if denotes -/V/(j^)(W^) in the case 
M; is GUE, 

EK] = np,S[y)) + O(^) and Var(r^) - + o(l)) logn. 

Hence, to establish Lemma [Sj it suffices to show that E[y„] — E\Y^ + o(l) and 
Var(y„) = Var(y^) + o(l). Below, we only deal with the variance, the argument 
for the expectation being similar and actually simpler. 

Set Ai — for i £ {1, . . . , n}. Notice that 

|Var(r„)-Var(i;:)| < ^ |(E[A,^,] -E[A,]E[A,]) - (Ei^^^] -E[A'JE[4])|. 

1< i,j<n 

Call an index i first class if E[Aj] > 1 — n^'^ or < n^'^ and second class otherwise. 

Note that if j is first class, then, for alH e {1. . . . , n}, 

\E[AA,]-E[AmAA\^0{n-^). 

Indeed, if E[Aj] < n^^, then both terms between the absolute value signs are 
0{n^^), so that |E[yliylj] - E[Ai]E[74j]| = 0{n^^). The other case can be brought 
back to this case by the identity 

|E[A,A,] - E[A,]E[A,]| = - E[B,]E[B,]|, 

where Bi is the complement of Ai. Consequently, 

(9) J2 linA^A,] E[A.]E[A,]) - (E[A:4] - E[A'JE[4])| = O(n-i). 

i or j 
1^^ class 

Theorem m shows that there are only O((logn)'-^'°siog") second class indices. 
Indeed, set rjn = ri~^e~^/^(logn)*-^'°siog" a^cj suppose first that i G {!,..., n} is 
such that t(i/n) < y ~ rjn'. 
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• if t{iln) > t{e), dH) is true for Wn- Then 

P(A, G /„) < P(|A, - tii/n)\ > 77„) < 

• if t{i/n) < t{e), choose j such that t{e) < t{j/n) < y — r]n (take e smaU 
enough and n large enough such that there is such a j). Then < Xj and 
P(A, £ In) = P(A, >y)< P(Aj >y)= P(Aj e /„) < 

Similarly one can show that if i G {1, . . . , n} is such that t{i/n) > y + r/n, then i is 
first class. 

As a consequence of this discussion, i G {1, . . . ,n} can only be second class if 
y — < t{i/n) < y + r/„. We need to count these possible i's. By definition of 
t{i/n), i = ^ fj^^"^ V^^dx. Thus, 



27r 



ry-v,^ ,j ry+v^ 

J-2 " " 2tt i_2 



1-2 ^" J-2 

In this case i belongs to an interval of length 
y+rin 2 



n 

27r Jy-ri.. 7r£ 



V4 - x^dx < __(logn)'^'°siog«_ 



Therefore, there are at most ^^(logn)'^i°si°8" + l = 0((log n)'^i°si°s") second 
class i's. 

Next, by Proposition [31 it is easily seen that if both i,j are second class, then 



n 



for some positive constant c. Since the number of such pairs is 0((logn)^'^'°s'°s")^ 

we have 

(10) 

J2 |(E[A^,]-E[A,]E[^,.])-(E[A^^;.]-E[^^]E[4.])l = 0(n-=(logn)2C^i°8'°sn). 

■i and j 
2'"* class 

To conclude, 

|Var(r„) - Var(i;:)| < ^ |(E[A,A,] - E[A,]E[A,]) - {E[A'^A'^] E[^^]E[4])| 

i or j" 
1=* class 



i and j 
2^'^ class 



so that jH) and ^ lead to 

|Var(y„) - Var(r,0| < 0(n"') + O(n-"(logn)2ci°siogn) ^ ^^^^^ 



as claimed. This shows that Var(y„) — 1^2^ + ''(l)j logn. As mentioned earlier, 
it may be shown similarly that E[y„] — npsc{I{y)) + o(l) and the proof of Lemma 
[S]is thus complete. □ 
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4. About real Wigner matrices 



In this section, we briefly indicate how the preceding resuhs for Hermitian ran- 
dom matrices may be stated similarly for real Wigner symmetric matrices. To this 
task, we follow the same scheme of proof, relying in particular on the corollary of 
Tao and Vu Four Moment Theorem (Proposition [3]) which also holds in the real 
case (cf. [6j). 

Real Wigner matrices are random symmetric matrices M„ of size n such that, for 
i < j, {Mn)ij are iid, with mean and variance 1, {Mn)u are iid with mean and 
variance 2. As in the complex case, an important example of real Wigner matrices 
is the case where the entries are Gaussian, giving rise to the so-called Gaussian 
Orthogonal Ensemble (GOE). 

The main issue is actually to establish first the conclusions for the GOE. This 
has been suitably developed by O'Rourke in j6] by means of interlacing formulas 



Theorem 6 (Forrester- Rains). The following relation holds between matrix ensem- 
bles: 



This statement can be interpreted in the following way. Take two independent 
matrices from the GOE, one of size n and the other of size n + 1. If we surimperpose 
the 2n+l eigenvalues on the real line and then take the n even ones, they have the 
same distribution as the eigenvalues of a n x n matrix from the GUE. 

Let / be an interval in M. Let be a GOE matrix and be a GUE matrix. 
W* and are the corresponding normalized matrices. The preceding interlacing 
formula leads to 



Relying on this result and on the GUE case, O'Rourke proved the following 
theorem: 



Theorem 7. Let be a GOE matrix. Set = ^M^. Then, for any y e 
(—2,2) and I{y) = [y, +oo), setting := Nit^y){W^), we have 



(cf. m). 



GUE„ = even(GOE„ U GOE„+i). 



. nNi{wn^nNi{w^)] + o{i) 

• \a.v{Ni{W^)) = 2Var(7V/(W^)) + 0(1), if Var(iV/(Wf )) c5o. 



E[y«] -n/,,c(/(y))+0(l) and 



Var(y«) = (-+o(l))logn, 




n- 



AA(0,1). 



oo 



Following exactly the same scheme as for complex Wigner matrices leads to the 
same conclusion: Theorem [7] is true for Wigner symmetric matrices, provided their 
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entries match the corresponding entries of GOE up to order 4 and satisfy condition 
(CO). 

The CLT for the eigenvalue counting function has been investigated as well for 
families of covariance matrices. The main conclusion of this work holds similarly in 
this case conditioned however on the validity of the Erdos-Yau-Yin rigidity theorem 
for covariance matrices. There is strong indication that the current approach by 
Erdos, Yin and Yau for Wigner matrices will indeed yield such a result. All the 
other ingredients of the proof are besides available. Indeed, Su (cf. [8j) carried 
out computations for Gaussian covariance matrices and proved both the CLT and 
the correct asymptotics for mean and variance. Tao and Vu in extended their 
Four Moment Theorem to such matrices. Arguing then as for Wigner matrices, we 
could reach in the same way a CLT with numerics for suitable families of covariance 
matrices. 
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